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COMPUTER-AIDED  DESIGN  OF  MILLIMETER  WAVE 


E-PLANE  FILTERS 


ABSTRACT 


This  study  presents  a  computer-aided  design  algorithm 
developed  for  a  class  of  E-plane  waveguide  filters.  The  filter 
structure  is  composed  of  an  integrated  finline  circuit  in  a 
waveguide.  The  circuit  consists  of  several  resonators  seperated 
by  inductive  strips.  The  analysis  portion  of  the  algorithm  is 
based  on  the  residue-calculus  technique  and  a  generalized  scat¬ 
tering  parameter  method.  It  is  mathematically  exact  and  numer¬ 
ically  very  efficient  because  the  convergence  is  guaranteed. 
Filters  designed  with  this  method  are  fabricated  and  tested 
at  Ka-band.  Good  agreement  with  design  is  obtained. 
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I .  INTRODUCTION 

Waveguide  filters  using  conventional  inductive  elements 
such  as  rods,  transverse  strips,  and  transverse  diaphrams  are  diffi¬ 
cult  to  produce  at  low  cost  and  large  quantities  because  of  their 
complicated  structures.  Microstrip  circuits  have  been  used  to  solve 
this  problem.  However,  these  circuits  are  typically  lossy,  espe¬ 
cially  at  millimeter-wave  frequencies.  A  number  of  low-loss  wave¬ 
guide  filters  have  been  proposed  recently  [1,2, 3, 4, 5]  that  are  com¬ 
posed  of  printed  circuits  inserted  in  the  waveguide  parallel  to  the 
E-plane  as  shown  in  Fig.  1.  The  circuits  are  designed  either  on 
a  metal  sheet  [1,2]  or  on  a  bilateral  finline  [3,4,5].  In  these 
structures,  the  slot  widths  are  equal  to  the  waveguide  height;  in 
effect,  the  circuits  consist  of  several  resonators  separated  by 
printed  Inductive  strips.  Since  the  structures  involve  only  straight 
line  shapes  and  are  amenable  to  photolithographic  techniques,  they 
are  highly  suitable  for  mass  production. 

In  the  previous  works,  most  of  the  design  procedures  are 
based  on  various  network-synthesis  techniques  where  the  inductive 
strips  are  represented  by  equivalent  T-networks.  The  equivalent 
circuit  parameters  were  obtained  either  by  measurements  [3]  or  by 
analytical  methods  such  as  mode  matching  techniques  [2,4]  and  varia¬ 
tional  methods  [1].  Although  these  procedures  are  accurate  enough 
for  some  designs,  the  equivalent  circuit  approach  neglects  the 
effect  of  higher-order  mode  couplings  between  the  strips.  The 


neglect  of  this  effect  causes  a  bandwidth  shrinkage  and  higher 
passband  ripples  in  the  filter  response.  To  alleviate  this  problem, 
Arndt,  et  al.  [5]  have  proposed  a  computer-aided  design  (CAD)  pro¬ 
cedure  that  takes  into  account  the  higher-order  mode  propagation 
for  design  of  the  bilateral  filters.  Their  analysis  is  based  on  a 
mode  matching  technique  that  involves  a  number  of  matrix  inversions. 
The  matrix  size  must  be  made  large  to  obtain  accurate  results,  and 
as  a  result,  the  analysis  becomes  inefficient  and  may  suffer  the 
problem  of  relative  convergence  [6]. 

In  this  study,  a  new  efficient  CAD  algorithm  is  developed 
for  a  class  of  E-plane  bandpass  filters  [7].  A  unified  treatment 
is  introduced  that  applies  not  only  to  purely  metallic  E-plane  cir¬ 
cuits  and  bilateral  finline  circuits  but  also  to  unilateral  and 
insulated  finline  circuits  (Fig.  2).  Since  the  analysis  portion 
of  this  algorithm  is  based  on  the  residue-calculus  technique  [8] 
and  a  generalized  scattering  parameter  method  [9,10],  it  is  mathe¬ 
matically  exact  and  numerically  very  efficient  because  the  conver¬ 
gence  is  guaranteed. 

In  the  next  chapter,  the  analysis  is  performed  on  a  loss¬ 
less  structure  with  infinitesimally  thin  conductive  sheets.  The 
residue-calculus  technique  is  first  applied  to  obtain  closed-form 
expressions  for  the  scattering  parameters  of  a  single  junction  cre¬ 
ated  by  a  semi-infinite  septum  in  the  waveguide.  Then  the  composite 
scattering  parameters  of  the  filter  structure  are  obtained  using 
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the  generalized  scattering  matrix  method  that  takes  into  account  the 
interaction  between  junctions  by  the  fundamental  mode  and  all  the 
higher-order  evanescent  modes.  Finally,  based  on  this  analysis  pro¬ 
cedure,  an  optimization  computer  program  is  used  to  design  filters 
for  a  given  performance  parameters.  Filters  designed  with  this 
method  are  fabricated  and  tested  in  Ka  band  [7] .  Good  agreement 
with  design  is  obtained. 

In  the  subsequent  chapter,  the  effect  due  to  the  thickness 
of  the  strips  is  discussed.  Based  on  the  results  in  the  previous 
analysis,  the  generalized  scattering  matrix  method  is  applied  again 
to  obtain  the  scattering  parameters  for  junctions  created  by  septa 
of  finite  thickness.  By  taking  into  account  the  metal  thickness, 
this  analysis  yields  results  in  good  agreement  with  experiments. 

II.  COMPUTER-AIDED  DESIGN  PROCEDURE 

The  CAD  program  consists  of  an  analysis  routine  and  an 
optimization  routine.  The  analysis  is  carried  out  in  three  steps. 
First,  a  limiting  case  of  a  single  junction  created  by  a  semi¬ 
infinite  septum  in  a  waveguide  is  analyzed.  This  structure  may  be 
analyzed  exactly  using  the  residue-calculus  technique  to  obtain  a 
closed-form  expression  for  the  scattering  matrix.  The  second  step 
is  to  calculate  the  scattering  parameters  for  a  finite-length  sep¬ 
tum.  This  is  done  by  placing  two  junctions  back-to-back  and 
utilizing  the  concept  of  the  generalized  scattering  matrix  that 


takes  into  account  the  interaction  between  junctions  by  the  propa¬ 
gating  modes  and  all  the  evanescent  modes.  Finally,  several  septa 
are  combined  to  form  a  filter  circuit.  The  total  scattering  para¬ 
meters  of  the  composite  structure  are  obtained  by  repeating  the  pro¬ 
cedure  in  the  last  step. 

Based  on  the  analysis,  an  objective  function  is  defined 
and  used  in  an  optimization  routine.  The  optimization  routine 
varies  the  input  parameters  until  the  objective,  i.e.,  the  desired 
value  of  the  insertion  loss  for  a  given  bandwidth,  is  achieved. 

Since  the  coupling  effects  and  frequency-dependent  properties  of 
the  dominant  mode,  as  well  as  higher-order  modes,  are  considered  in 
the  analysis,  the  design  procedure  yields  accurate  results.  In 
addition,  a  wider  range  of  filter  performance  may  be  obtained 
through  the  optimization  procedure. 

A.  Step  1  —  Semi-Infinite  Septum 

Let  us  consider  an  infinitesimally  thin  semi-infinite  sep¬ 
tum  in  a  waveguide  lnhomogeneously  filled  with  four  dielectric  slabs 
as  shown  in  Fig.  3.  This  is  a  generalized  structure  for  which 
values  of  the  thickness  and  dielectric  constant  of  each  layer  can 
be  chosen  to  represent  the  structures  in  Fig.  1  and  2.  (Note  that 
for  the  case  of  bilateral  finlines,  a  magnetic  wall  is  Introduced  at 
the  center  of  the  waveguide  and  only  the  region  above  the  magnetic 
wall  is  considered.  Such  a  region  represents  the  upper  half  of  the 
bilateral  finline  due  to  the  symmetry  of  the  structure.)  It  is  then 
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subdivided  into  three  regions  at  the  junction  z  ■  0,  with  each 
region  consisting  of  a  uniform  waveguide  which  possesses  known  eigen- 
modes  and  corresponding  propagation  constants. 


When  a  TE^-type  wave  strikes  on  the  junction,  the  excited 


fields  are  composed  of  the  TE^-type  waves  since  there  are  no  struc¬ 


tural  variations  in  the  y-direction.  Therefore,  the  existing  modal 


fields  consist  of  only  E^,  and  components.  The  total  fields 


in  each  region  are  then  expanded  in  terms  of  the  orthonormal  eigen- 


modes  as 
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e<„  -  •  lp  + i  *<„♦.„<*) 


iy  ip  ip 


in  in 
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(  1) 


where  {A.  }  are  unknown  coefficients  and  A'  is  the  amplitude  of 


the  incident  wave  from  the  i-th  region;  the  expressions  for  the 


mode  function  d>in (x)  and  the  propagation  constant  y^  are  given  in 


Appendix  A  for  different  structures.  The  continuity  condition  at 


z  -  0  requires  that  the  tangential  electric  field  (E  )  and  magnetic 

1  3E  y 

field  (H  ■  - - )  be  continuous.  This  leads  to  the  following 
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relations: 
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From  these  relations  a  set  of  equations  involving  the  unknown  coef¬ 
ficients  is  derived  by  making  use  of  the  mode  orthogonality  property 
This  is  done  by  multiplying  (  2)  by  ant*  and  *nte“ 

grating  with  respect  to  x  from  o  to  b  and  from  b  to  a,  respectively. 
The  results  are: 
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This  Infinite  set  of  equations  is  to  be  solved  simultaneously  for 

the  unknown  coefficients  {A.  }. 

in 

An  alternative  set  of  equations  may  also  be  derived  by 
multiplying  the  equations  (  2)  by  (x)  and  integrating  with  res 


pect  to  x  from  o  to  a.  The  resulting  equations  are: 


|>  Wtjfea  fail 


‘  *-*  V  ^  **  > #' 


a;  F  “  A.  F  A!  G  »  A,  G  A 

2r  p _  +  y  2p  n  +  _ 3£_£ _  +  c  _ ?n  n  „  -  _m 

-Y2p  Ylm  n*l  Y2n  ^lm  -^3p  Ylm  n*l  Y3n  Ylm  Hm 


(  5b) 


Thus  far,  the  mode-matching  formulation  has  led  to  two 
sets  of  equations,  (  4)  and  (  5) .  Since  the  existence  of  the 

layered  dielectrics  does  not  affect  the  edge  condition  of  the  sep¬ 
tum,  the  equations  are  of  the  same  form  as  that  for  a  bifurcated 
waveguide  described  in  [8]  and  the  scattering  parameters  of  the 
junction  are  obtainable  by  the  residue  calculus  technique. 

If  there  is  only  one  incident  wave  coming  from  region  I, 

i.e.,  Al  *  Al  ■  0,  (  4)  becomes 
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Consider  the  integrals 
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where  the  contour  C  is  an  infinitely  large  circle  in  the  complex  w- 
plane,  traversed  in  the  counterclockwise  direction.  Let  the  func¬ 
tion  f(w)  satisfy  the  following  conditions: 


1.  f(w)  is  an  analytic  function  of  w  except  for  simple 
poles  at  w  ■  y^n»  n  -  1,2,3,...  and  at  w  ■  -y^  . 

2.  f(w)  has  simple  zeroes  at  w  *  -y„  ,  and  w  =  -y_  , 

jm 

xq  *  *  • 

3.  The  residue  of  f(w)  at  v  *  -Y-  is  equal  to  A*  H  . 

ip  *  iP  P 

4.  f(w)  has  algebraic  behavior  at  infinity.  Specifically, 
f(w)  -  0(w~3^2). 

In  view  of  condition  4,  the  integrands  in  (  7)  are 

-5/2 

0(w  )  on  the  contour  C;  consequently,  these  integrals  are  equal 

to  zero.  Expressing  the  integrals  in  terms  of  series  of  residues 
on  f(w)  and  equating  the  results  to  zero  we  have 
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for  m  *  1,2,3,...,  where  the  condition 


has  been  applied.  A  comparison 
diately  shows  that. 


3  of  f(w),  Resf(-Y^)  = 
of  (1-8)  and  (1-6)  inme- 
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The  unknowns  (A.  }  are  readily  determined  if  f(w)  is  known, 
in 

At  this  stage  the  reason  for  enforcing  the  four  conditions 
on  f(w)  becomes  obvious.  The  pole  structure  given  in  condition  1 
was  chosen  to  make  the  residue  series  Involving  {y  }  identical  to 
the  left  hand  side  of  (  6),  while  the  assigned  locations  of  the 

zeroes  of  f(w)  in  condition  2  ensure  that  no  contribution  comes  from 
the  poles  at  w  *  ~Y2m  an<*  ~Y3m  contour  integrals.  Condition 

3  brings  the  dependence  of  the  incident  wave  into  the  final  solu¬ 
tion  and  may,  therefore,  be  regarded  as  a  normalization  of  f (w) . 
Condition  4  not  only  ensures  that  the  contour  integrals  are  zero, 
but  also  allows  unique  determination  of  the  solution  by  means  of 
the  specific  large-argument  asymptotic  behavior  of  f(w).  Condition 

4  is  related  to  the  edge  condition  [11].  This  point  will  become 
clear  when  f(w)  is  to  be  constructed  for  solutions. 

The  construction  of  f(w)  is  the  key  step  in  the  residue- 
calculus  method.  First,  the  zero  and  pole  conditions  in  1  and  2 
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In  a  similar  manner,  assumption  of  a  single  incident  wave  A2p^2p^ 
or  A3p^3p^x^  1°  equation  (  4)  and  repetition  of  the  procedure  yield 

the  expressions  for  S^fapp)  or  S^(m,p),  respectively.  They  are 
summarized  in  the  following: 
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Notice  that  in  the  preceding  equation  (  14)  the  expressions 


(Y,  -  Y,  ) ,  (y~  -  Yo  ) ,  and  (y  -  y,  )  are  introduced  to  cancel 

lm  im  zp  Zp  op  jp 

the  corresponding  poles  in  the  infinite  products.  The  notation 
S^(m,p)  stands  for  the  scattering  coefficient  of  the  m-th  order 
mode  in  region  i  due  to  the  p-th  order  mode  excitation  incident 
from  region  j.  This  set  of  equations  is  thus  valid  for  any  mode. 

An  alternative  set  of  solutions  for  S^(m,p)  can  be 
obtained  by  solving  the  set  of  equations  in  (  5) .  The  results  pro¬ 

vide  a  valuable  checking  tool  in  numerical  calculations.  They  are, 
therefore,  summarized  in  Appendix  C. 


B.  Step  2  —  Finite-Length  Septum 

In  step  1,  we  have  characterized  the  single  junction  by 
means  of  a  three-port  scattering  matrix  that  has  nine  elements. 

Each  of  the  elements  is  also  a  matrix  of  infinite  dimensions  that 
corresponds  to  the  infinite  number  of  eigenmodes.  With  the  know¬ 
ledge  of  the  scattering  parameters  for  a  single  junction,  the 
generalized  scattering  matrix  technique  is  applied  to  obtain  a  two- 
port  scattering  matrix  as  shown  in  Fig.  4.  Suppose  now  that  the 
TE-type  wave  from  region  I  is  incident  upon  junction  A.  At  this 
junction,  fields  are  partly  reflected  back  into  region  I  and  partly 
transmitted  into  regions  II  and  III.  After  traveling  a  distance  d, 
part  of  the  wave  transmitted  into  regions  II  and  III  is  reflected 
back  and  part  is  transmitted  into  region  I*  at  junction  B.  This 


process  continues  until  the  intensity  of  the  reflected  wave  dies 
out.  The  multiple-reflection  phenomenon  between  junctions  A  and  B 
is  implied  in  a  matrix  manipulation  which  yields  the  scattering 
matrix  for  a  finite  length  septum. 

In  Fig.  4,  and  represent  the  scattering  matrices 
for  isolated  junctions  A  and  B  (i.e.,  for  semi- inf inite  s ep turns ) , 
respectively.  From  equation  (  14)  we  obtain  the  value  of  the 

elements  in  both  matrices.  Since  the  characteristics  of  these  two 
junctions  are  essentially  the  same  except  for  the  opposite  orienta¬ 
tion,  if  we  define  the  port  numbers  of  the  networks  as  shown  in  Fig.  4, 
we  have 


S,  =  S  =S 
b  a 

where  S  is  the  scattering  matrix  defined  by  (14). 
We  now  define  the  transmission  matrix 
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where  [I]  and  [OJ  are  the  identity  matrix  and  the  zero  matrix, 

respectively;  [T  ]  and  [T„]  are  diagonal  matrices  of  infinite  size 

-V21d/2 

whose  diagonal  elements  are  *  e  and  T3ii  =  e  for 


i  from  1  to  00 ,  respectively.  In  effect,  [Tg]  and  [T^]  represent 
the  wave  propagating  (for  propagation  modes)  or  attenuating  (for 
evanescent  modes)  for  a  distance  of  d/2  in  guided  regions  II  and  III, 
respectively.  The  combination  of  the  transmission  matrix  along  with 
S  and  S.  results  in  S'  and  S.' : 
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A  further  combination  of  S'  and  S'  then  yields  the  composite  scat- 
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where 


s . .  =  l  S',  S' 
~ij  kt2  lk  kJ 


i,j  *  1,2,3 


that  represents  the  septum  of  length  d,  viewed  from  I  at  z  =  0  and 
from  I'  at  z  *  d. 


C.  Step  3  —  Filter  Structure 

For  the  final  step  of  the  analysis,  several  finite  length 
septa  are  cascaded  with  appropriate  separations  to  form  a  filter 
circuit.  This  is  done  by  using  the  same  techniques  used  in  step  2 
described  above.  The  procedure  depicted  in  Fig.  5  shows  the  cas- 

<\j  %  % 

cading  of  two  septa  represented  by  S  and  S,  .  The  comination  of  S 

3  U  S 

and  the  transmission  matrix  T,  results  in  matrix  S'  which  is  fur- 

i  a 

ther  combined  with  to  yield  S,  the  total  scattering  matrix  of 
two  septa  with  a  separation  of  £. 
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Derivation  of  S  parameters  for  cascading  septa. 
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This  process  is  repeated  to  obtain  the  scattering  parameters  of  a 
filter  structure  consisting  of  any  number  of  septa. 

D.  Numerical  Considerations 

Based  on  the  above  analysis,  a  computer  program  is  devel¬ 
oped  to  calculate  the  total  scattering  parameters  for  a  filter  struc 
ture  with  given  design  parameters.  This  analysis  program  is  to  be 
used  later  in  an  optimization  routine.  Since  the  scattering  para¬ 
meters  of  the  semi- inf inite  septum  depend  only  on  the  waveguide 
size  and  dielectric  properties  that  are  constant  factors  in  a 
design,  the  calculations  in  step  1  need  to  be  done  only  once  at  each 
frequency  point.  These  calculations  may  be  regarded  as  pre¬ 
calculations  to  the  optimization  routine;  therefore,  the  optimiza¬ 
tion  procedure  requires  only  the  calculations  in  steps  2  and  3. 

Although  the  analysis  is  formally  exact,  in  practice  the 
infinite  products  in  (  14)  and  the  infinite  dimensions  of  the  ele¬ 

ments  in  the  scattering  matrices  need  to  be  truncated.  To  guarantee 
the  reliability  of  numerical  results,  the  degree  of  convergence  in 
the  numerical  calculations  is  studied.  The  numerical  results  show 
that  the  infinite  products  in  (  14)  converge  very  fast.  Table  1 


gives  two  typical  calculations  of  S^(l,l)  and  ^^(l,  1)  as  the 
number  of  terms  in  the  product  is  increased.  Although  fairly 
accurate  results  may  be  obtained  by  using  the  first  few  terms  (<  10), 
at  least  300  terms  are  required  to  achieve  an  accuracy  of  0.1  per¬ 
cent.  The  calculation  of  the  products  here  is  efficient  since  only 
readily  known  propagation  constants  are  involved. 

For  most  applications,  the  fundamental  mode  is  below  cut¬ 
off  in  the  narrow  waveguide  section  (regions  11  and  111)  and  all  of 
the  higher  order  modes  are  evanescent  in  the  wide  waveguide  (region 
1) .  Therefore,  only  the  first  few  accessible  modes  [12]  are  essen¬ 
tially  responsible  for  the  coupling  between  resonators.  For  most 
cases,  only  two  or  three  modes  are  required  to  describe  accurate 
scattering  behavior.  Typical  examples  of  S^(.l,l)  for  septa  of 
lengths  d  ■■  1  mm  and  5  mm,  calculated  using  up  to  5  modes,  are  shown 
in  Table  2.  It  is  clear  that  more  modes  need  to  be  used  for  a 
short  septum  at  a  higher  frequency.  Notice  that  the  number  of 
modes  considered  here  indicates  the  dimensions  of  the  matrices  to 
be  manipulated  in  steps  2  and  3.  The  calculation  includes  a  number 
of  complex  matrix  multiplications  and  inversions.  This  analysis 
is  efficient  mainly  due  to  the  fact  that  the  matrix  size  is  kept 


small. 
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E.  Optimization 

The  analysis  algorithm  is  then  used  by  an  optimization 
computer  program  [13]  that  varies  the  input  parameters  until  a 
desired  value  of  the  Insertion  loss  for  a  given  bandwidth  is  obtained. 
This  is  done  by  defining  an  error  (objective)  function  F(x)  as  fol¬ 
lows: 


F(x) 


(  20) 


where  Ng  and  are  the  number  of  frequency  sampling  points  in  the 
stopband  and  passband,  respectively;  Lg  and  are  the  minimum  stop- 
band  and  maximum  passband  attenuation  levels,  respectively;  | 

is  the  calculated  value  of  the  filter  attenuation.  For  a  given 
thickness  of  the  dielectric,  the  parameters  x  to  be  optimized  are 
the  septum  lengths  and  resonator  lengths. 

To  prevent  the  program  from  approaching  non-realizable  or 
non-physical  results,  upper  and  lower  bounds  of  each  parameter  are 
specified 
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The  constrained  parameters  are  then  reduced  into  essentially  uncon¬ 
strained  ones  by  the  transformation  [14] 
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where  -«®  <  <  00  but  only  solutions  within  the  range 

0  <  cot_1x^  <  n  (  23) 

are  allowed.  This  transformation  has  a  penalizing  effect  upon  the 
parameters  in  the  vicinity  of  the  upper  and  lower  bounds. 

For  the  optimization  procedure,  the  first  two  modes  in 
each  subregion  are  considered;  the  final  design  values  are  calcu¬ 
lated  with  three  modes.  The  total  computing  time  for  optimizing  a 
five  resonator  filter  was  about  2  minutes  on  a  CDC  Dual  Cyber  170/ 
750. 

F.  Experiments 

In  an  earlier  publication  [15],  it  has  been  pointed  out 
that  the  Q-factor  in  insulated  finlines  is  inferior  to  that  of 
bilateral  finlines.  Therefore,  the  design  of  filters  using  insul¬ 
ated  finlines  will  not  be  considered  in  this  study.  Further,  the 
E-plane  circuits  usually  employ  a  relatively  thick  metal  sheet; 
these  circuits  will  be  discussed  in  the  next  chapter. 

In  the  present  study,  bandpass  filters  in  bilateral  and 
unilateral  finline  configurations  are  designed  and  fabricated  on 
RT/Duroid  substrate.  Some  typical  designs  for  Ka-band  applications 
are  shown  in  Figures  6  and  8.  Notice  that  the  resonators  are 
very  weakly  coupled  by  wide  septa.  In  such  cases,  the  septum  can¬ 
not  be  represented  by  simple  shunt-inductance  equivalent  circuits. 
Instead,  it  has  to  be  represented  by  an  equivalent  T  network  [1], 
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Since  the  equivalent  circuits  are  useful  in  conventional  circuit 
design  and  are  helpful  for  good  initial  guesses  in  the  optimization 
procedure,  this  topic  is  discussed  briefly  in  Appendix  D. 

The  five-resonator  filter  design  in  Fig.  6  is  a  bilateral 
structure  in  which  the  hatched  portion  indicates  metal  on  both  sides 
of  the  substrate.  The  center  frequency  is  38.85  GHz;  the  0.1  dB 
ripple  bandwidth  is  1.1  GHz;  and  the  skirt  selectivity  is  -40  dB 
at  500  MHz  away  from  either  edge  of  the  passband.  Measured  charac¬ 
teristics  of  this  filter  are  shown  in  Figure  7.  Computed  responses 
are  also  plotted  for  comparison.  Shapes  of  the  calculated  and  mea¬ 
sured  response  curves  agree  extremely  well.  Measured  insertion 
loss  in  the  passband  is  less  than  1.3  dB,  of  which  0.3  dB  is  asso¬ 
ciated  with  the  test  fixture.  The  measured  center  frequency  is 
shifted  to  the  lower  side  by  about  170  MHz,  which  is  less  than  0.5Z 
of  the  center  frequency. 

Fig.  8  shows  a  photograph  of  a  three-resonator  bilateral 
finline  filter  printed  on  a  10-mil  thick  RT/Duroid  substrate.  The 
design  specifications  for  this  fitler  were  the  same  as  those  for  the 
last  one  except  for  a  decreased  skirt  selectivity  of  -20  dB.  The 
measured  results  in  Fig.  9  show  that  the  insertion  loss  is  less 
than  0.8  dB  in  the  passband,  including  the  fixture  loss.  The  cal¬ 
culated  and  measured  frequency  responses  again  agree  well;  however, 
the  measured  center  frequency  is  shifted  down  by  about  300  MHz. 


Three-Resonator  Bilateral  Finline  Filter 
Resonator  Length: 

£^  =  £^  =  2.6559  mm 
£ 2  ~  2.6187  mm 

Septum  Width: 

d,  =  d,  =  1.7230  mm 
1  4 

&2  “  d^  *  5.3628  mm 


8.  Three-resonator  Ka-band  bilateral  filter 


Bandpass  filters  have  also  been  developed  in  unilateral 


finline  configurations.  A  typical  design  for  a  three-resonator 
filter  has  the  following  specifications  and  dimensions: 

Resonator  lengths:  =  £3  =  2.1568  mm;  =  2.0879  mm 

Septum  widths:  d^  =  d^  =  2.3325  mm;  d^  =  d^  =  6.8718  mm 

Passband:  39.0  ^  39.7  GHz,  Chebyshev  response  with  0.1 
dB  ripple 

Skirt  selectivity:  -20  dB  at  500  MHz  away  from  passband 

edge 

The  measured  results  in  Fig.  10  show  about  2  dB  insertion  loss  in 
the  passband  and  about  200  MHz  downward  shift  for  the  center  fre¬ 
quency. 

In  the  experiments  performed,  the  measured  center  fre¬ 
quencies  are  always  lower  than  the  calculated  values.  The  devia¬ 
tions  range  from  170  MHz  to  300  MHz.  One  cause  of  this  shift  is 
that  the  filter  structure  assembled  is  different  from  the  one 
analyzed.  The  assembled  structure  has  10-mil-deep  grooves  machined 
into  the  top  and  bottom  walls  to  support  the  substrate  in  the  wave¬ 
guide.  An  experimental  study  on  this  effect  shows  that  within  the 
quarter-wavelength  range,  the  deeper  the  groove  is  made,  the  more 
the  center  frequency  will  shift  down.  Another  possible  source  of 
the  discrepancy  is  the  finite  thickness  of  the  metal  fin.  The 
thickness  effect  is  minor  at  this  frequency  since  the  thickness  of 


the  copper  cladding  on  the  substrate  is  about  0.5  mil.  However, 
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10.  Frequency  response  of  a  3-resonator  unilateral  filter 
(a)  Insertion  loss  characteristics,  (b)  Return  loss 
characteristics . 


it  will  become  an  important  factor  for  higher  frequency  applications, 
such  as  W-band  circuits.  These  effects  are  discussed  in  the  next 
chapter. 

From  the  experiments,  it  has  been  observed  that  the  inser¬ 
tion  loss  of  a  unilateral  filter  is  always  greater  than  that  of  a 
bilateral  filter  for  similar  design  specifications.  One  possible 
explanation  is  that  the  dielectric  loss  of  the  unilateral  filter  is 
greater,  while  the  conductor  loss  is  about  the  same  in  both  cases. 

In  the  case  of  bilateral  finlines,  little  energy  can  penetrate  into 
the  very  narrow  space  between  the  strips.  Therefore,  the  energy  is 
mainly  coupled  through  the  wider  empty  waveguide  section.  In  the 
unilateral  case,  the  energy  is  coupled  about  evenly  through  the 
empty  guide  and  the  partially  dielectric-filled  guide  region. 
Obviously,  more  dielectric  loss  is  encountered  in  this  case.  By 
the  same  reasoning  we  can  surmise  that  the  conductor  loss  is  about 
the  same  for  both  cases.  Although  the  bilateral  structure  has  cop¬ 
per  cladding  on  both  sides  of  the  substrate,  there  is  essentially 
no  current  flowing  on  the  inner  surface  facing  the  dielectric.  It 
is  thus  equivalent  to  a  conducting  strip  with  current  flowing  on 
the  outer  surface  of  both  sides  as  in  the  unilateral  case.  Of 
course,  these  statements  are  strictly  conjectural;  further  study  is 
needed  to  verify  them.  This  conjecture  would  also  be  applicable  to 
the  case  of  insulated  finlines  versus  bilateral  finlines. 


Another  point  worth  mentioning  is  the  sensitivity  of  the 
filter  structure  to  mechanical  tolerances.  Calculations  show  that 
for  both  cases,  the  center  frequency  shifts  about  100  MHz  if  the 
fabrication  error  in  resonator  length  is  1  mil.  Also,  if  the  dielec 
trie  substrate  is  shifted  in  the  broadwall  direction  (x-direction) 
from  its  centered  position,  the  filter's  center  frequency  falls. 

This  effect  is  much  more  noticeable  with  the  unilateral  structure. 
Therefore,  the  bilateral  structure  shows  more  promise  from  a  con¬ 
structional  viewpoint. 

III.  FILTERS  WITH  SEPTA  OF  FINITE  THICKNESS 

In  the  previous  chapter,  we  have  introduced  a  design 
procedure  for  a  class  of  E-plane  waveguide  filters.  The  metal 
septa  used  in  the  circuits,  either  stand-alone  or  printed  on  a 
substrate,  were  considered  to  be  infinitesimally  thin.  This  assump¬ 
tion  is  applicable  only  when  the  thickness  of  the  septa  is  small 
compared  to  the  wavelength,  typically  less  than  0.3%.  When  the 
purely  metallic  sheet  is  used,  or  when  the  operating  frequency  is 
high,  this  criterion  is  difficult  to  satisfy,  and  the  thickness 
effect  must  be  considered  in  the  analysis  for  an  accurate  design. 

In  this  chapter,  the  design  procedure  is  modified  to  take 
into  account  the  thickness  of  the  metal  sheets.  To  this  end,  step 
1  of  the  previous  analysis  is  replaced  by  the  analysis  of  a  semi¬ 
infinite  septum  of  finite  thickness.  The  rest  of  the  design  pro¬ 
cedure  is  kept  unchanged. 


A.  Finite-Thickness  Septum 

Fig.  11a  shows  the  generalized  structure  to  he  analyzed, 
similar  to  that  in  Fig.  3  except  for  the  septum  of  thickness  t. 
Again  the  excitation  field  is  assumed  to  be  of  TE  type.  Despite 
the  modal  expansions  for  the  fields  on  either  side  of  the  interface 
z  =  0,  matching  of  the  transverse  components  does  not  lead  to  a  set 
of  equations  of  the  form  solvable  hy  the  residue-calculus  technique. 
However,  the  modal  nature  of  the  fields  in  the  three  regions  sug¬ 
gests  description  of  the  junction  in  terms  of  generalized  scat¬ 
tering  matrices.  The  scattering  matrices  may  be  derived  via  the 
multiple-reflection  method  for  a  suitable  auxiliary  geometry. 

The  auxiliary  geometry  appropriate  for  the  thick-septum 
junction  is  that  of  Fig.  lib.  The  conducting  wall  due  to  the 
thick  septum  is  recessed  into  region  III  and  two  additional  regions 
IV  and  V  are  created  for  convenience  of  analysis.  The  original 
structure  can  be  recovered  by  reducing  both  6^  and  6^  to  zero. 

Notice  that  in  this  new  structure,  the  characteristics  of  junctions 
at  z  =>  0  (junction  A)  and  z  =  6^  (junction  B)  are  essentially  the 
same  as  that  of  the  thin-septum  junction  whose  exact  scattering 
matrix  description  was  given  in  the  previous  chapter.  Let  the  scat¬ 
tering  matrices  for  the  Junction  A  have  the  subscript  a,  i.e., 

(i,j  «  1,2,4);  and  let  those  for  the  junction  B  have  the  subscript 
b,  i.e.,  (i,j  *  4,5,3).  The  scattering  matrices  for  the  com¬ 

posite  junction, (i,j*  1,2,3),  will  have  no  subscript.  Con¬ 
sideration  of  the  multiple-reflection  phenomena  and  reduction  of 
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Fig.  11.  Semi-infinite  septum  of  finite  thickness.  (a)  Geometry 
(b)  Auxiliary  structure. 


6  ,  §2  to  zero  yields  the  following  representations  for  the  compo¬ 
site  scattering  matrices: 


511  =  Sall  +  Sal4fl  '  Sb44Sa44J  Sb44Sa41 

512  =  Sal2  +  Sal4fI  "  Sb44Sa44J  Sb44Sa42 

513  =  Sal4[l  ~  Sb44Sa44^  Sb43 

521  =  Sa21  +  Sa24[l  "  Sb44Sa44^  Sb44Sa41 

522  =  Sa22  +  Sa24[I  ~  Sb44Sa44^  Sb44Sa42 

523  *  S24fI  “  Sb44Sa441  Sb43 

531  "  Sb34[I  "  Sa44Sb44]  Sa41 

532  =  Sb34tI"  Sa44Sb44]  Sa42 

533  =  Sb34l1  "  Sa44Sb44]  Sa44Sb43 
where 


(  24a) 

(  24b) 

(  24c) 
(  24d) 

(  24e) 

(  24  f) 

(  24g) 

(  24h) 

(  241) 


H 


B.  Numerical  Considerations 

The  computation  of  equation  (  24)  involves  a  number  of 

matrix  inversions  and  multiplications.  The  computation  time  can 
be  significantly  reduced  by  truncating  the  scattering  matrices  at 
the  smallest  size  that  yields  suitably  accurate  results.  A  conver¬ 
gence  study  shows  that  for  most  of  the  cases  only  3x3  matrices 
are  required.  Some  typical  examples  are  shown  in  Tables  3  and  4 
where  the  dominant  scattering  parameters  of  the  junction  created  by 
a  finite-thickness  septum  are  calculated  using  matrices  of  sizes 
ranging  from  1  x  1  to  6  x  6.  It  is  clear  that  in  this  case  a  2  x  2 
matrix  calculation  yields  results  accurate  to  the  third  digit. 

C.  Thickness  Effect 

From  a  physical  point  of  view,  the  larger  the  scattering 
area  of  an  obstacle,  the  more  energy  is  reflected.  However,  in  a 
bifurcated  waveguide,  the  change  in  the  scattering  area  (i.e.,  the 
thickness  of  the  septum),  does  not  affect  the  amplitude  of  the 
reflection  coefficient  since  all  of  the  energy  is  reflected  no 
matter  how  thin  the  septum  is.  It  is  the  phase  angle  of  the  reflec¬ 
tion  coefficient  that  is  affected  by  the  change  in  thickness.  The 
computed  phase  angle  for  a  septum  located  at  the  center  of  the  wave¬ 
guide  is  shown  in  Fig.  12  as  a  function  of  the  normalized  thick¬ 

ness  t/a.  It  shows  that  the  phase  angle  increases  with  an  increasing 
thickness  of  the  septum.  If  we  replace  the  bifurcated  waveguide 
section  by  an  equivalent  empty  waveguide  of  length  £  short  circuited 


at  the  very  end,  the  thickness  effect  means  that  i  decreases  with 
an  increasing  t.  This  is  because  when  t  increases,  the  narrow- 
waveguide  sections  become  smaller  and  less  energy  is  stored  in  them. 

The  effect  of  thickness  on  the  scattering  characteristics 
of  a  finite-length  septum  is  more  explicit.  It  is  no  longer  true 
in  this  case  that  all  the  electromagnetic  energy  is  reflected; 
instead,  part  of  the  energy  is  transmitted  through  the  cut-off  wave¬ 
guide  by  way  of  evanescent  mode  coupling.  It  is  obvious  that  the 
increase  of  septum  thickness  would  cause  more  reflection  and  less 
transmission.  This  phenomenon  is  well  shown  by  the  results  sum¬ 
marized  in  Table  5. 

Finally,  the  thickness  effect  of  the  septum  on  the  per¬ 
formance  of  the  filter  circuits  is  demonstrated  by  a  filter  designed 
on  a  purely  metallic  sheet.  The  performance  of  the  bandpass  filter 
is  first  optimized  on  a  5-mil  thick  sheet.  Then  using  the  same 
design  parameters,  the  filter  response  is  calculated  for  the  filter 
circuits  by  varying  the  thickness  of  the  metal  sheet.  The  results 
are  plotted  in  Fig.  13.  The  most  noticeable  effect  is  that  tue 
center  frequency  shifts  upward  approximately  1  MHz  per  1-mil 
increase  in  thickness.  This  effect  is  mainly  due  to  the  fact  that 
the  increase  of  septum  thickness  effectively  shortens  the  equiva¬ 
lent  length  of  the  resonators  formed  between  two  septa,  and  thus 
causes  the  resonant  frequency  to  shift  upward.  In  addition  to  this 
effect,  there  are  some  other  effects.  For  example,  the  passband 
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Fig.  13. 


Thickness  effects  on  filter  performance  (second  filter  design 
in  Table  6.) 


ripple  increases  perceptibly  and  the  bandwidth  shrinks  slightly 


due  to  changes  in  coupling  between  resonators. 


D.  Design  Examples 

Several  E-plane  filters  employing  septa  with  finite  thick¬ 
ness  have  been  designed  and  tested  at  microwave  frequencies  by 
Konishi,  et  al.  [1]  and  Tajima,  et  al.  [2].  The  validity  of  this 
analysis  is  checked  by  using  their  design  data  to  calculate  the 
filter  responses.  The  results  are  shown  in  Fig.  14  in  comparison 
with  their  measurements.  Good  agreement  is  observed. 

Filters  for  Ka-band  and  W-band  applications  are  designed 
with  pure  metal  structures.  Typical  design  parameters  are  tabulated 
in  Table  6.  There  are  certain  advantages  inherent  to  these  struc¬ 
tures:  First,  there  is  no  dielectric  loss  involved.  Second,  the 
grooves  for  supporting  the  circuits  in  the  waveguide  are  totally 
filled  with  metal;  therefore,  the  error  due  to  the  grooves  is 
reduced.  It  is  expected  that  these  filters  are  superior  to  those 
made  in  unilateral  and  bilateral  finline  configurations. 


IV.  CONCLUSIONS 

A  class  of  low-loss  waveguide  filters  can  be  designed 
using  the  unified  CAD  algorithm  developed  in  this  study.  The  filter 
structure  is  composed  of  a  printed  circuit  inserted  in  a  waveguide 
parallel  to  the  E-plane.  This  circuit  is  designed  either  on  a 
purely  metallic  sheet  or  on  various  finline  structures.  It  consists 


INSERTION  LOSS  l dB)  ~  INSERTION  LOSS  ( <JB  ) 


Frequency  response  of  E-plane  filters  using  metal  sheets 
of  finite  thickness. 


of  several  resonators  separated  by  inductive  strips;  such  a  con¬ 
figuration  is  highly  suitable  for  low-cost  mass  production. 

The  analysis  portion  of  the  CAD  algorithm  is  based  on  the 
combination  of  two  methods,  viz,  the  residue-calculus  technique  and 
the  generalized  scattering  matrix  method.  The  former  is  first 
applied  to  obtain  closed-form  expressions  for  the  scattering  para¬ 
meters  of  a  junction  created  by  a  semi-infinite  septum  in  the  wave¬ 
guide.  Based  on  this  Information,  the  latter  is  then  applied  to 
obtain  the  composite  scattering  parameters  for  a  filter  structure. 
This  analysis  is  mathematically  exact  and  numerically  very  effi¬ 
cient  because  the  convergence  is  guaranteed. 

Filters  have  been  designed  and  tested  in  Ka  band  using 
bilateral  and  unilateral  configurations.  Measured  characteristics 
agree  well  with  theoretical  data.  Extensive  experiments  show  that 
the  bilateral  constructure  is  more  promising  than  the  unilateral 
structure. 

The  thickness  effect  of  the  metal  sheet  is  also  dis¬ 
cussed.  Neglecting  this  effect  in  the  design  causes  a  shift  in  the 
center  frequency  and  degrades  the  performance  of  the  filter.  The 
analysis  taking  into  account  the  metal  thickness  has  been  incor¬ 
porated  in  the  design  procedure,  which  is  believed  to  be  useful  for 
design  of  filters  above  Ka  band. 
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APPENDIX  A.  EXPRESSIONS  FOR  EIGENFUNCTIONS 

The  function  has  to  satisfy  the  second-order  ordinary 


differential  equation 


32<j> 


3x 


in  ,2 

*—  -  k  <p.  =  0 

2  xn  in 


(25) 


and  the  boundary  conditions  at  x  =  0  and  x  =  a.  The  wave  number  in 
the  x  direction,  kxn,  is  related  to  the  propagation  constant  Y^  via 
the  dispersion  relation 


2  J  u2 

Y  -  k  =  e  k 
n  xn  r  o 


(26) 


where  kQ  is  the  free-space  wave  number.  The  expressions  for  $  and 
of  different  geometries  are  given  in  the  following  : 
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,  0  <  x  <  b 


b  <  x  <  a 


APPENDIX  B.  EXPRESSIONS  FOR  H  ,  F  ,  AND  G 

n  n’  n 


The  coefficients  H  ,  F  ,  and  G  are  related  to  the  structure 
n  n*  n 


geometry  and  the  field  distribution  in  regions  I,  II,  and  III 


respectively.  They  are  derived  from  the  overlapping  integrals 


/  *lm(x)  *2n(x)  dx  "  ~2~  2~ 

°  Y2n  -  *lm 


(5  2  a) 


I  *lm(x)  ^3n(x)  dx  =  T~  n  T 
b  Y3n  -  Ylm 


(3  2  b) 


where  <()^n(x)  ,  ^2n^X^  *  and  ^3n^X^  are  tbe  n-tb  order  orthonormal 


eigenfunctions  for  the  E^  fields  in  regions  I,  II,  and  III,  res¬ 


pectively  (refer  to  Fig.  3).  The  expressions  for  F^,  and 


of  different  geometries  are  given  in  the  following. 
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Bilateral  Finline  Structure 
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where  B  and  E  are  defined  in  equations  (28-2  9) 
n  n 
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C.  ALTERNATIVE  SET  OF  SCATTERING  PARAMETERS 


-YlpHm  eL<Ylm+  Ylp)  (Ylp-  Ylp) 


Y,  H 
lm  p 


<YiP+ 


8<Ylm,  Yi.) 


lm’  'lp' 


(36a  ) 


^l£  L(Y2m+  Ylp)  (Ylp"  YlP)(Y2m"  Y2m)  ,  w 

HF  (Yl_+  Y,J  8(Y2m’YlP)  (36b) 


p  m 


rlp  ’  2m 


eL(Y3.+  Ylp>  <Ilp-  Y1P)<Y3„-  V3m) 

HG  (Yl^+  Y-J„)  8lY3m’Ylp'  {'3bc) 


Li£ 


p  m 


'lp'  T3m' 


-H  ] 
m 

*p  L(Y2p+  Yln) 

2 

/  \ 

(  36d ) 

2y: 

c. 

Lm 

<Y2p+Ylm> 

-F 

_ 2. 

L(Y2m+ 

V 

(Y2m‘ 

Y2m) 

g(Y2m’Y2p) 

(  36e) 

F 

m 

e 

^Y2m+ 

Y2p> 

-F 
_ 2 

L(Y3m+ 

Y2p} 

<Y3m~ 

Y3m} 

8(Y3m’Y2p) 

(  36f) 

G 

m 

e 

(Y3m+ 

V 

-H»Cp  „L<Y3p+  V 

1 

(  36g) 

2Ylm 

( V  1,i») 

8(Ylm’Y3p) 

-G 
_ 2 

L(Y2m+ 

Y3p} 

<Y2m~ 

Y2mJ 

g(Y2m’Y3p) 

(  36h) 

F 

m 

e 

(Y2m+ 

Y3p} 

-G 
_ 2 

L(Y3m+ 

V 

(Y3m‘ 

Y3m} 

g(Y3m’Y3p) 

(  36i) 

G 

m 

e 

<Y3m+ 

Y3p} 

n  (Ym+  0t)(Y2n+  B)CY3n4  B) 

n-1  (Yln-  S)(Y2n-  «)(Y3n-  a) 


g(  a,B  ) 


APPENDIX  D.  EQUIVALENT  CIRCUIT  OF  A  SEPTUM 


The  equivalent  T-network  for  a  finite-length  septum  is  shown 

in  Fig.  15  .  The  values  of  the  normalized  reactances  x  and  x  are 

s  p 

to  be  calculated  from  the  scattering  parameters  of  a  septum  of 
length,  say,  D.  Let  the  scattering  matrix  of  the  fundamental  mode 
be 


(37  ) 


For  this  reciprocal  lossless  two-port  device,  the  following  relations 
hold: 


S12  “  S21 
S11  =  S22 

|SUI2  +  is12l2  -  1 


These  relations  provide  useful  error-checking  criteria  for  the 
numerical  calculations. 

If  we  convert  the  scattering  matrix  into  impedance  matrix  and 
equate  it  to  the  impedance  matrix  of  the  equivalent  T  network,  we 
have 


[Z]  “  [I  +  S]  [I  -  S] 
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(39) 


P  s  -  p 

where  I  is  the  identity  matrix.  By  equating  the  elements  in  the 
matrix  on  both  sides,  we  obtain  x  and  x  as: 
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